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Abstract. The generation of parallel electric fields by the propagation of ion cyclotron waves 
(with frequency 0.3 u>d) in the plasma with a transverse density inhomogeneity was studied. 
Using two- fluid, cold plasma linearised equations, it was shown for the first time that, in this 
particular context, En generation can be understood by an analytic equation that couples E\\ 
to the transverse electric field of the driving ion cyclotron wave. It was proven that the minimal 
model required to reproduce the previous kinetic simulation results of En generation [Tsiklauri 
et al 2005, Genot et al 2004] is the two-fluid, cold plasma approximation in the linear regime. 
By considering the numerical solutions it was also shown that the cause of En generation is 
the electron and ion flow separation induced by the transverse density inhomogeneity. We also 
investigate how En generation is affected by the mass ratio and found that amplitude attained 
by En decreases linearly as inverse of the mass ratio mi/m e . For realistic mass ratio of rrii/m e = 
1836, such empirical scaling law, within a time corresponding to 3 periods of the driving ion 
cyclotron wave, is producing En = 14 Vm -1 for solar coronal parameters. Increase in mass ratio 
does not have any effect on final parallel (magnetic field aligned) speed attained by electrons. 
However, parallel ion velocity decreases linearly with inverse of the mass ratio mi/m e . These 
results can be interpreted as following: (i) ion dynamics plays no role in the En generation; (ii) 
En cc 1/m.i scaling is caused by the fact that w<j = 0.3uj C i oc 1/rrii is decreasing with the increase 
of ion mass, and hence the electron fluid can effectively "short-circuit" (recombine with) the 
slowly oscillating ions, hence producing smaller En. 
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The generation of parallel electric fields in inhomogeneous plasmas is a generic topic, 
which is of interest in a variety of plasma phenomena such as particle acceleration in Solar 
and stellar flares Fletcher(2005), auroral acceleration region and current sheets in the 
Earth magnetosphere (see refs. in Song & Lysak(2006)), laboratory plasma reconnection 
experiments Yamada & et al(1997) and many more. In situ and remote observations of 
accelerated particles often show parallel electric fields in localised double layers, charge 
holes or U-shaped voltage drops. 

The issue of Eu generation by the propagation of ion cyclotron waves in the plasma 
with a transverse density inhomogeneity is discussed in more detail in Tsiklauri(2007). 



We start from two-fluid, cold (ignoring thermal pressure) plasma linearised equations 
Krall & Trivelpiece(1973): 
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d t E = cV x B - 4irne(Vi - V e ). (0.4) 

Hereafter subscripts under d denote partial derivative with respect to that subscript. 
Uniform, background magnetic field, B is in z-dircction. Density profile is specified as a 
ramp, n(x) = uq (l + 3exp [-[(a; - 100<5)/(20<5)] 6 ] ) in which the central region (along in- 
direction, i.e. across z), is smoothly enhanced by a factor of 4, and there are the strongest 
density gradients having a width of about 205 around the points x = 815 and x = 1195. 
Here 5 — c/uj pe is the (electron) skin depth, which is a unit of grid in our numerical 
simulation. We use 2.5D description meaning that we keep all three, x, y, z components 
of all vectors, however spatial derivatives d/dy = 0. 

In order to derive the equation that describes Eu = E z generation, we write Eqs.(l)-(4) 
in x, y, z component form. Omitting details of the calculation we present the final result: 

(d 2 tt - c 2 d 2 xx + w 2 pi + cj 2 pe ) £|| = -c 2 d 2 zx E x . (0.5) 

Also, a similar calculation enables us to obtain the equation describing the dynamics of 
driving transverse electric field E x on an ion cyclotron wave: 

(d 2 t -c 2 d 2 zz +w 2 pi +u 2 pe )E x = 

-c 2 d 2 zx E\\ - ujl^mi/ejuJdViy - uo 2 pe {m e / e)uo ce V ey . (0.6) 

Note that Eq.(6) also describes the feedback of the generated En on the driving transverse 
electric field E x (see the first term on the right-hand-side). Here the notation is standard: 
ujpe — ^/47rn e 2 /m e and tu P i = \J '47m e 2 jvm are electron and ion plasma frequencies; 
w c(e.i) — e -So/( m ( e ,i) c ) are respective cyclotron frequencies. 

In order to solve Eqs.(l)-(4) numerically we use the following normalisation: t = 
tu> pi }, V XtVtZ = V x . VtZ c, E x , VtZ = E x>y ^ z (m e cuj pe /e) — E Xj y iZ Eo, B x , ViZ = B x ^_ z Bq, and 
(x,y,z) = cj uj pe {x , y , z) . In what follows we omit tilde on the dimensionlcss quantities. 
The (x, z) simulation 2D box size is 200(5 x 2500(5. Since we fix background plasma num- 
ber density at 10 9 cm" 3 (typical value for the solar corona), lo p& is then 1.784 x 10 9 rad 
s" 1 and the simulation box size is 33.6 m in x- and 420.5 m in z-dircction. Bq was fixed 
at 101.5 Gauss (typical value for the solar corona), which gives uj ce /oj pe — 1. mi/m e ratio 
was varied as: 45.9, 91.8, 183.6 to 262.286 (realistic one is 1836). These values correspond 
to 1/40, 1/20, 1/10 and 1/7-th of the realistic value respectively. This yields respectively: 
io ci /uj pl = Ba/icy/lirrnrrn) = V A /c = l/^m,/m e = 0.148, 0.104, 0.074 and 0.062 for 
x ^ 70 and x ^ 130 (realistic u> C i/oj p i = Va/c is 0.023). Here parameters are similar 
to e.g. Tsiklauri & et al(2005a), Tsiklauri & et al(2005b), except for far more realis- 
tic mass ratios. Note that the simulation parameters are still somewhat artificial. Full 
kinetic, Particle-In-Cell (PIC) simulations employed in Tsiklauri & et al(2005a), Tsik- 
lauri & et al(2005b) or in gyro-kinetic approach which uses guiding centre approximation 
for electrons, whilst retaining ion particle-like dynamics Gcnot & et al(1999), Genot & 
et al(2004) are computationally challenging. Thus, in those studies rather modest mass 
ratios e.g. 16 were used. Note also, that since here we do not need to resolve electron 
thermal motions as we are only studying electromagnetic part of the problem {E\\ gen- 
eration) our unit of spatial grid size is 5 = c/Lu pei the (electron) skin depth. While in 
full kinetic, PIC simulation Tsiklauri & et al(2005a), Tsiklauri & et al(2005b) the unit 
of grid has to be A = v t h,e/uj pe . Since in a PIC simulation typically vth,e/c = 0.1, 
in the present, two-fluid approach an equivalent to PIC numerical simulation requires 
(5/ A) 2 = (c/v t h,e) 2 = 10 2 less grid points, thus it can be 100 times faster. For compari- 
son a single run for mass ration 16 in Tsiklauri & et al(2005a), Tsiklauri & et al(2005b) 
takes about 8 days on parallel, 32 dual-core 2.4 GHz Xeon processors, similar run with 
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Figure 1. Contour plot of E z = Eu at time t = 16000/o; ce . Here m,i/m e = 262. 



mass ratio of 262 would have taken 4 months. The numerical run presented here for the 
mass ratio of 262 takes 4 days with only one processor. 

We solve relativistic version of Eqs.(l)-(4) numerically with a specially developed and 
tested FORTRAN 90 code which uses 4-th order centred spatial derivatives and 4-th 
order Runge-Kutta time marching. Although Alfven speeds considered are at most « 15 
% of the speed of light for rrii/m e — 45.9, relativistic effects were included. The simplest 
option becomes available in the linear regime. In Krall & Trivelpiece(1973), appendix 
I, paragraph 5, it was shown that the relativistic equation of motion of a particle with 
charge q and the rest mass too can be written as 

dt 7TO0 

where 7 = (1 - V 2 /c 2 )-^ 2 . As can be seen from the latter equation, in the linear 
regime, it coincides with either Eq.(l) or (2) after substituting m e ^ — > Jennie,!, where 
7e,j = (1 — Ve i/c 2 )" 1 ^ 2 - Naturally, such simplified approach is only valid when there are 
no flows in the unperturbed state Vq = 0. As can be seen below, largest attained velocities 
in the simulation are those of electrons, and these do not exceed 3 % of speed of light. 
Thus, relativistic corrections play only a minor role. It should be noted, however we still 
retain the displacement current in Eq.(4). Note, also that the gradients in the code are 
resolved numerically to an appropriate precision (20 grid points across each gradient.) 

Initially all perturbations are set to zero, and we start driving the 2 = 1 cell with the 
transverse magnetic fields of the form B y = — 0.05 sin^i) (l-O — exp[— (t/(3.125w~ 1 )) 2 ]) 
and B x = -0.05cos(w d t) (1.0 - exp[-(t/(3.125w~ 1 )) 2 ]) . As in Tsiklauri & et al(2005a), 
Tsiklauri & et al(2005b), we fixed tod at 0.3 u> C i (to avoid ion-cyclotron damping play- 
ing any role). (l.O — exp[— (t/(3.125w^ 1 )) 2 ]) factor ensures that these driving B± fields 
ramp up to their maximal values in time t — 3.125cj ( ^ 1 . Such driving with B± of 5% of 
the background Bo excites circularly polarised ion-cyclotron (IC) waves, these waves are 
often misquoted as Alfven waves Genot & et al(1999), Genot & et al(2004), Tsiklauri & 
et al(2005a), Tsiklauri & et al(2005b). Although, in the frequency range oj <§; io ci both 
left and right polarised IC waves tend to an Alfven wave branch in the dispersion rela- 
tion Krall & Trivelpiece(1973), at frequencies u> ~ 0.3w C i the correct term ion-cyclotron 
wave instead should be used. In the considered problem E x and B y are both components 
of Alfven (ion cyclotron) wave, so these can be used interchangeably. Conventionally, 
Alfvenic and IC waves are more associated with magnetic field oscillation. It is in ki- 
netic, Particle-In-Cell simulations driving by electric field is required because particles 
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Figure 2. Contour plot of (V e z ~ Viz) oc j z at time t = 16000/o; ce . Here m,i/m e — 262. 

respond to electric, rather than magnetic fields. In the two-fluid simulation this is not a 
requirement. 

The generated at the left edge {z = 1) IC waves propagate both in the directions 
of positive and negative z's. However, because of the periodic boundary conditions used 
(applied on all physical quantities) IC wave that travels to the direction of negative z's (to 
the left) re-appears on the right edge of the simulation box. The driving B± dynamics 
is not presented here (see for details Tsiklauri(2007)), but it resembles closely to the 
previous kinetic simulation results Genot & et al(1999), Genot & et al(2004), Tsiklauri & 
et al(2005a), Tsiklauri & et al(2005b). As in all previous phase- mixing simulations Alfven 
velocity is a function of the transverse (to the background magnetic field) coordinate, 
x, i.e. Va = Va{x) oc 1/ yjn{x). Thus the IC wave middle portion travels slower than 
the parts close to the simulation box edge. This creates progressively strong transverse 
gradients and hence smaller spatial scales. If resistive effects are included (these are 
absent here) , such a configuration usually produces greatly enhanced dissipation and IC 
wave amplitude decays in space as oc exp(— z 3 ) Tsiklauri & et al(2005a), Tsiklauri & 
et al(2005b). 

The E\\ = E z field snapshot at t = 16000/o; ce is shown in Fig.(l). 

is generated only in the regions of density gradients i.e. around x = 81 and x = 119 
lines. This can be explained by analysing right-hand-side (RHS) of Eq.(5). En = at 
t = everywhere, however it can only be generated in the regions where d x E x ^ 0. The 
latter is true only in the density gradient regions where E x becomes progressively oblique 
propagating. Thus, Eq.(5), derived here for the first time, correctly explains the simulated 
process of En generation by IC waves. Also, note that E\\ amplitude at time t = lBOOOw^ 1 
attains value of 3 x 10~ 5 . This is somewhat smaller value than the one obtained in the full 
kinetic (PIC) simulation Tsiklauri & et al(2005a), Tsiklauri & et al(2005b). This is due to 
the different mass ratios: in the the kinetic (PIC) simulation mi/m e — 16, but here it is 
262. In dimensional units this E\\ corresponds to about 0.003 statvolt cm -1 or 90 V m -1 , 
i.e. in such electric field electrons would be accelerated to the energy of ~ 10 keV over 
the distance of 100 m. Note, however, that the generated E\\ is oscillatory in space and 
time. The typical values of the Dreicer electric field on the corona is few x 10~ 3 V m _1 
Tsiklauri(2006b), which implies the obtained Eu in our model exceeds the Dreicer value 
by at least four orders of magnitude, perhaps enabling the electron run away regime. 
This would imply that our model is more relevant to the acceleration of solar wind, 
rather than solving coronal heating problem. Essentially acceleration of electrons would 
dominate over the heating as such. However, this seems uncertain because electron and 
ion fluid separation cannot go on forever, and some sort of discharge should eventually 
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take place. At any rate, similar kinetic simulations have shown Genot & et al(2004) (see 
their Fig. (11)) that that wave energy is converted into particle energy on timescales of 
lO 3 ^^ 1 (mind that the latter number is likely to be dependent on the mass ratio mj/m e ). 

In Fig. (2) we present (V ez — V iz ) which is proportional to j z , the parallel (electron and 
ion flow separation) current. It can be seen from this figure that (V ez — Vi Z ) attains mod- 
erate values of w 0.03c. Genot & et al(2004) stated the importance of charge separation 
before. However the cause of E\\ generation is actually electron and ion flow separation 
(see below). The latter is quite different from the electrostatic effect of charge separation, 
which is inherently a plasma kinetic effect. Electron and ion flow separation is a fluid- like 
(non-kinetic) effect, because their distribution functions remain Maxwellian at all times. 

In Fig. (3) we present of the final attained parallel electric field amplitude (within 
3 periods of the driving ion cyclotron wave) as a function of mass ratio, i.e. E* = 
max.(\E z (x, z,tfi na i)\) vs. ■m i jvn e . We see that the amplitude attained by E\\ decreases 
linearly with inverse of the mass ratio mj/m e . The ai-range in Fig. (3) is mj/m e = 30 — 
1836, so that rightmost point of the dashed line enables us to grasp E\\ for the case of 
realistic mass ratio (i.e. 1836). We thus gather that E\\ = 0.0085/1836 = 4.630 x 10~ 6 
which is 4.630 x 10~ 6 x E Q = 4.7 x 10~ 4 statvolt cm" 1 or 14 Vm" 1 . Note that it is likely 
that the actual value will be even smaller if the dissipativc effects (resistivity, viscosity) 
will be taken into account. However, these are know to be very small in the solar coronal 
plasma. Thus, we do not expect a big change in this result. 

In Fig. (4) we show what trends in generation of parallel (magnetic field aligned) 
electron and ion flow as a function of the mass ratio. It follows from Fig. (4) that 
an increase in mass ratio does not have any effect on final parallel speed attained 
by electrons. However, the parallel ion velocity decreases linearly with inverse of the 
mass ratio m,/m e . The ratio of final attained ion and electron flow amplitudes (within 
3 periods of the driving ion cyclotron wave) as a function of mass ratio, i.e. V* = 
max(\V iz (x,z,tf ina i)\)/ max(|V ez (x, z, t fina i)\), shows clear scaling of l/(m,/m e ). 

The conclusions that follow from the collective analysis of Figs.(3)-(4) initially may 
seem counterintuitive. On one hand maximal attained E\\ amplitudes drop off as l/rrij 
(Fig. (3)). On the other hand, electron flow maximal attained amplitudes do not depend 
on mi (they all are circa 0.03c, see Fig. (4)), while ion flow maximal attained ampli- 
tudes are much smaller (0.0001c— 0.00065c) than that of electrons and drop off as l/m* 
(Fig. (4)). Thus one might expect that more massive ions should produce a bigger Ej| 
(since separation of electron and ion fluids is the source of E\\ and that separation is 
expected to be largest in the case of more massive ions, as they are slower). In fact, this 
is what would be expected if the polarisation drift produced by the driving IC wave is the 
cause of parallel electric field generation Genot & et al(1999), Genot & et al(2004). The 
latter two references use the following polarisation drift current, j± = (m^i/ B 2 )dE±/dt. 
This equation implies that E» then should increase with oc m;. However, in Fig. (3) we 
see completely opposite E\\ oc 1/m^ scaling. These results can be interpreted (reconciled) 
as following: (i) ion dynamics plays no role in the E\\ generation, i.e. polarisation drift 
has no effect in contrary to the claims of Genot & et al(1999), Genot & et al(2004); (ii) 
decrease in the generated parallel electric field amplitude with the increase of the mass 
ratio mi/m e is caused by the fact that aid = 0.3w C i cx 1/m, is decreasing, and hence the 
electron fluid can effectively " short-circuit" (recombine with) the slowly oscillating ions, 
hence producing smaller E\\ which also scales exactly as 1/mj. 

It should be noted that since here we use two-fluid approach the generated E\\ can- 
not change the distribution function, which obviously remains Maxwellian, while in the 
previous kinetic simulation of a similar system it produced bumps in the distribution 
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Figure 3. A log-log plot of the final attained parallel electric field amplitude (generated within 
3 periods of the driving ion cyclotron wave) as a function of mass ratio. Data points correspond 
to E* = max(\E z (x,z,t fina i)\) for m;/m e = 45.9, 91.8, 183.6 and 262.286. The dashed line 
corresponds to the fit 0.0085/(m;/m e ). 
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Figure 4. A log-log plot of the final attained electron and ion flow amplitudes (generated 
within 3 periods of the driving ion cyclotron wave) as a function of mass ratio. Solid squares 
correspond to max(| V ez (x, z, tf ina i)\), open circles to max(\Vi z (x, z, tf ina i)\), while open squares 
show their ratio, i.e. V* = max(|l / iz (x, z, tfi na i)\)/ max(\V ez (x, z, t/» n oi)D- The solid line show a 
fit 1/x — l/(rrii/m e ). 



function as the electrons residing on the magnetic field lines with the density gradients 
get efficiently accelerated (see e.g. Fig. (4) in Tsiklauri & et al(2005b)). 

The generation of En is a generic feature of plasmas with the transverse density inho- 
mogcncity and in a different context this was known for decades in the laboratory plasmas 
Cross & Miljak(1993), Ross & et al(1982). Also, it should be emphasised that the two 
fluid description in the context of parallel electric field generation has been used before 
Goertz & Boswell(1979). Moreover, Lysak(1990) has explored similar equation. It can be 
easily shown that our Eq.(5) can be reduced to Eq.(44) from Lysak(1990), by neglect- 
ing the displacement current and writing it in Fourier-transformed form. At first sight, 
this may seem to diminish the importance of our result. However, first, the inclusion of 
the displacement current mathematically means introduction of the time-dynamics (df t 
term in Eq.(5)) which is crucial for the correct description of E\\ generation by IC wave 
driving in the transversely inhomogeneous plasma; second, generally the importance of 
the inclusion of the displacement current in this context has been recently stressed Song 
& Lysak(2006); and third, because of transverse inhomogeneity of plasma it is impos- 
sible to use Fourier-transform in the transverse spatial coordinate. It should be noted 
when plasma density is homogeneous no Eu generation takes place, in our model; and 
this is corroborated both by numerical simulations (not presented here) and agrees with 
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the Eq.(5) (when n — const, the RHS of Eq.(5) is zero at all times as E x does not 
propagate obliquely). However, the relation of this claim to the observations is not en- 
tirely straightforward. Observations in the Earth's magnetosphere have for many years 
noted that perpendicular gradients of the perpendicular electric field are commonly seen 
during auroral zone crossings Mozer & et al(1977). These are sometimes seen in the 
presence of perpendicular (transverse) density gradients; and the plasma density typi- 
cally decreases above parallel potential drops Carlson & et al(1998). The question here 
is what is generated first E\\ or E± (since obviously both are coupled via Eq.(5))? Our 
model also correctly reproduces the previous kinetic results Genot & et al(1999), Genot 
& et al(2004), Tsiklauri & et al(2005a), Tsiklauri & et al(2005b) that only electrons are 
accelerated (along the background magnetic field), while ions practically show no accel- 
eration. This applies only on time scales considered i.e. t ^ 16000/o; ce . However, in the 
Earth's auroral zone a significant amount of ion acceleration takes place Horwitz(1982) 
and Alfvenic Poynting flux is directly correlated with ion outflow in the auroral zone 
Strangeway & et al(2000). Based on this one can conjecture that either processes re- 
ported by these observations acted over time scale t » 16000/kJ ce so that ions due to 
their large inertia have enough time to accelerate or a mechanism other then proposed 
here is responsible for the ion acceleration. It should be noted that because of oscillatory 
nature of the obtained .Eii , it can possibly act as yet another mechanism for interpreting 
the peculiar hard x-ray (> 25 keV) solar flare, which is believed to be produced by a 
non-thermal electron beam Ofman & Sui(2006), Nakariakov & et al(2006). 
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